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^p^l Abstract 

2 . The X-ray Imaging Spectrometer (XIS) on board the Suzaku satellite is an X-ray CCD camera system 

' that has superior performance such as a low background, high quantum efficiency, and good energy res- 

Ci . olution in the 0.2-12 keV band. Because of the radiation damage in orbit, however, the charge transfer 
inefficiency (CTI) has increased, and hence the energy scale and resolution of the XIS has been degraded 
since the launch of July 2005. The CCD has a charge injection structure, and the CTI of each column and 

^ ■ the pulse-height dependence of the CTI are precisely measured by a checker flag charge injection (CFCI) 

V/^ ' technique. Our precise CTI correction improved the energy resolution from 230 eV to 190 cV at 5.9 keV 

■ in December 2006. This paper reports the CTI measurements with the CFCI experiments in orbit. Using 
^ ■ the CFCI results, we have implemented the time-dependent energy scale and resolution to the Suzaku 

, calibration database. 

Q^ ' Key words: instrumentation: detectors — techniques: spectroscopic — X-ray CCDs 

o 

OO 

■ 1. Introduction The increase of the CTI is due to an increase in the num- 
ber of charge traps at defects in the lattice structure of 

. ^ , After the first successful space flight use of the X-ray silicon made by the radiation. Since the trap distribution 

' charge coupled device (CCD) of the SIS (Burke et al. 1993) is not uniform, it would be best if we could measure the 

JH ' on board ASCA, the CCD has been playing a major role CTI of each pixel as Chandra ACIS (Grant et al. 2004). 

. - - ' in imaging spectroscopy in the field of X-ray astronomy. In the case of the XIS, however, it is impossible to mea- 

However, the charge transfer inefficiency (CTI) of X-ray sure the CTI values of all the pixels, mainly because the 

CCDs increases in orbit due to the radiation damage; the onboard calibration sources do not cover the entire field of 

CTI is defined as the fraction of electrons that are not sue- view of the XIS. Therefore, we use the CTI of each column 

cessfuUy moved from one CCD pixel to the next during the to correct the positional dependence of the energy scale, 
readout. Since the amount of charge loss depends on the The XIS is equipped with a charge injection struc- 

number of the transfers, the energy scale of X-ray CCDs ture (Prigozhin et al. 2004; Bautz et al. 2004; LaMarr 

depends on the location of an X-ray event. Furthermore, et al. 2004; Prigozhin et al. 2008) which can inject an ar- 

there is a fluctuation in the amount of the lost charge, bitrary amount of charge in arbitrary positions. Using 

Therefore, without any correction, the energy resolution this capability, we can precisely measure the CTI of each 

of X-ray CCDs in orbit gradually degrades. In the case column (Nakajima et al. 2008). By applying the column- 

of the X-ray Imaging Spectrometer (XIS) (Koyama et al. to-column CTI correction, the positional dependence of 

2007) on board the Suzaku satellite (Mitsuda et al. 2007) the CTI corrected energy scale is greatly reduced, and 

launched on July 10, 2005, the energy resolution in full the over-all energy resolution is also improved (Nakajima 

width at half maximum (FWHM) at 5.9 keV was -140 et al. 2008). 

eV in August 2005, but had degraded to ~230 eV in In Nakajima et al. (2008), the results of the CTI cor- 

December 2006. rection was mainly based on the ground-based charge in- 
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jection experiments. In-orbit measurements were limited 
within one year after the launch. This paper reports more 
comprehensive and extended in-orbit experiments up to 
two years after the launch. The results arc based on the 
data with the normal full window mode (Koyama et al. 
2007) without a spaced-row charge injection^ (Uchiyama 
et al. 2009), and have been implemented to the Suzaku 
calibration database and applied to all the data obtained 
with the same mode. All the errors arc at the Icr confi- 
dence level throughout this paper unless mentioned. 

2. X-Ray CCD with the Charge Injection 
Capability 

2.1. Checker Flag Charge Injection 

The XIS is the set of four X-ray CCD camera systems. 
Three sensors (XIS 0, 2, and 3) contain front-illuminated 
(FI) CCDs and the other (XIS 1) contains back illumi- 
nated (BI) CCD. The XIS 2 sensor became unusable on 
November 9, 2006. Therefore there are no data for XIS 2 
after that day. 

The detailed structure of the CCD has been provided 
in Koyama et al. (2007). In this paper, we define a "row" 
and a " column" as a CCD line along the ActX and ActY 
axes, respectively (see figure 3 in Koyama et al. 2007). In 
the imaging area, the ActX value runs to 1023 from the 
segment A to D, while the ActY value runs from to 1023 
from the readout node to the charge injection structure. 

The charge injection structure lies adjacent to the top 
row {ActY = 1023) in the imaging area. We can inject 
charges from ^ 50 e~ to ~ 4000 e~ per pixel; the equiva- 
lent X-ray energy ranges from ~0.2 keV to ~15 keV. 

A charge packet generated by an incident X-ray is trans- 
ferred to the readout node, then is converted to a pulse- 
height value. We define PHq to be the original pulse 
height generated by the X-ray. In the real case, the read- 
out pulse height of the packet [PH') is smaller than PHq, 
because some amount of charges is lost during the trans- 
fer. To measure the charge loss, we have to know both 
PHo and PH' . However, we can usually measure only 
PH' , and hence it is difficult to obtain PHq. 

Prigozhin et al. (2004) and Prigozhin et al. (2008) re- 
ported a technique to solve this problem by the charge 
injection method, and Nakajima et al. (2008) applied this 
technique to the XIS. We briefly repeat by referring fig- 
ure 3 in Nakajima et al. (2008). First, we inject a "test" 
charge packet into the top CCD row {ActY = 1023). 
Then, after the gap of a few rows, five continuous packets 
are injected with the same amount of charge of the test 
packet. The former four packets are called "sacrificial" 
charge packets, while the last one is called a "reference" 
charge packet. The test packet loses its charge by the 
charge traps. On the other hand, the reference packet 
does not suffer from the charge loss, because the traps are 
already filled by the preceding sacrificial packets. Thus we 
can measure the charge loss by comparing the pulse-height 

^ After October 2006, the spaced-row charge injection (Uchiyama 
et al. 2009) has been a normal observation mode, and hence we 
should use different correction method. 



values of the reference charge {PH^ef) and the test charge 
{PHtest)- The relation between sacrificial charge pack- 
ets and reference charge packets is described in Gendreau 
(1995). We can obtain a checker flag pattern by these in- 
jected packets in the X-ray image (right panel of figure 3 
in Nakajima et al. 2008) because of the onboard event- 
detection algorithm (Koyama et al. 2007). Therefore in 
this paper, we call this technique a "checker flag charge 
injection (CFCI)." 

2. 2. Formulation of the CTI 

A charge packet in the XIS loses its charge during 

(a) the fast transfer (24 /is pixel"^) along the ActY 
axis in the imaging area, (b) the fast transfer along the 
^ciFaxis in the frame-store region, (c) the slow transfer 
(6.7 ms pixel"^) along the ActY axis in the frame-store 
region, (d) the fast transfer to the readout node along the 
ActX axis. The CTI depends on many parameters such as 
the transfer speed and the number density of the charge 
traps (Hardy et al. 1998). The frame-store region is cov- 
ered by the shield and is not exposed to the radiation 
directly. Furthermore, the pixel size of the frame-store re- 
gion (21 /imxl3.5 ^m) is different from that of the imag- 
ing area (24 /xmx24 fim). Thus the number of traps per 
pixel may be different between the imaging area and the 
frame-store region. Then we assumed that the four trans- 
fers have different CTI values. We examined the transfer 
(d) by using the calibration source data taken in April 
2007, and found no significant decrease of the pulse height 
along the ActX axis. We, therefore, ignore the charge loss 
in the transfer (d). 

We define that i is the transfer number in the imaging 
area {i=ActY+l; here, ActY is a coordinate value where 
an incident X-ray generates a charge packet). Then the 
relation between PH' and PHo is expressed as, 

PH'{t) = PHq{1 - Cay{l - Ct)'°'^-\l - Ce)\ 

-Pi?o{l-«(Ca-Cfc + C,)- 1024c,,}, (1) 

where Cq, Cb and Cc are the CTI values in the transfers (a), 

(b) , and (c), respectively. Here we used the fact that the 
CTI values are much smaller than 1. Thus we can separate 
the charge loss into i-dependent component (the second 
term in the right-hand side of equation 1) and constant 
component (the third term). We therefore substitute the 
CTI with CTI 1 (the former component) and CTI 2 (the 
latter component), which have the CTI values of ci = 
Ca — Cb + Cc and C2 = c^, respectively. Then equation 1 can 
be written as 

PH'{i) ^ PHo{l-ici- 1024:C2). (2) 

Since the CTI values depend on the amount of trans- 
fer charge which is proportional to the pulse height, we 
assume the CTI is described by a power function of the 
pulse height (Prigozhin et al. 2004) and expressed as 

ci = fci(Pi7o)~'^ and C2 = k2{PHQ)-'^, (3) 

where fci and fc2 are scale factors for the CTI 1 and CTI 2, 
and the index f3 is common to the CTIl and CTI 2. 
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3. Measuring the CTI with the CFCI 

We have conducted the CFCI experiments six times in 
orbit. Effective exposure time for each experiment ranges 
from a few to ~20 ks. The equivalent X-ray energy of the 
injected charge packets ranges from ~0.3 keV to ~8 keV. 
Since June 2006, we injected various amounts of charge in 
one experiment. The log is summarized in table 1. 

In the CFCI experiments, the test charge is injected to 
the row at ActY=1023 (i = 1024), and hence Pil'(1024) = 
PHtcst- The reference charge should be equal to the origi- 
nal charge which does not suffer from the charge loss, and 
hence PHo = PH^cf- Then equation 2 can be written as 

^^ + ^^= 1024(,^-Pi^j- 

We determined ci 4- C2 by measuring the ratio 
PHtast/ PPti-cf for each column. From equation 3, we can 
obtain the relation in the CFCI experiments as, 

Cl+C2^iki+k2)PHj. (5) 

The index /? and ki + k2 were derived by fitting equation 5 
to the values of ci -I- C2 obtained with the CFCI experi- 
ments with multiple amount of charge injections (multi- 
ple PHicis; Log Number 3-6 in table 1). The mean and 
standard deviation of the best-fit /3 of equation 5 averaged 
over each sensor are shown in figure 1 . The mean value of 
(3 shows no time variation, and the time averaged values 
of XIS 0, 1, and 3 are 0.31, 0.22, 0.15 for XIS 0, 1, and 3, 
respectively. As for XIS 2, there was only one data point, 
and we obtained /3 = 0.34. 

If a charge packet has a volume proportional to the 
number of electrons and is spherically symmetric, the 
probability that one electron encounters a charge trap is 
proportional to the cross section of the charge packet. In 
this case, we can expect {PH-ccf — Pi/tcst) PH^lf ■ From 
equation 4, the CTI value is proportional to {PH^cf — 

— 1/3 

PHtest) / PHicii and hence oiPH^^^ . Thus the simple 
model is roughly consistent with the observed (3 values. 

Using the results of all the CFCI experiments (Log 
Number 1-6 in table 1) and above determined (3 values, 
we then re-estimated fci and k2 separately. In this process, 
we assumed ki/k2 is equal to Ci/c2 which was estimated 
by the 6.4 keV line from the Sgr C region to be 0.67 and 
1.5 for the FI and BI CCD, respectively (Nakajima et al. 
2008). From this fci, we can obtain the final value of ci. 

Figure 2 shows an example of the distribution of ci in 
July 2006. We can see significant column-to-column dis- 
persion. Figure 3 shows the change of ci from July 2006 
to September 2007. We can see that the CTI values of 
all columns increased, but the increasing rate was differ- 
ent from column to column. The results of figures 2 and 
3 indicate that the CTI correction at the column level 
is strongly required. In figure 4, we show the column- 
averaged Ci value as a function of time. 

Since the CFCI experiments were only sparsely con- 
ducted (see table 1), we interpolate the ci and C2 values 
for the observations of inter-CFCI epochs. As for the de- 
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termination of the CTI values before the first CFCI ex- 
periment, see appendix. 

4. Results of the CTI Correction and the Energy 
Calibration 

A CTI correction, which is the conversion of PH'{i) 
to PHo, is made with equation 2, where ci and C2 are 
calculated from equation 3 by using the fci, fe, and /3 
values determined in section 3. 

Emission Lines for the Calibration 

We used the emission lines from the onboard calibration 
sources, the Perseus cluster of galaxies, and the supernova 
remnant IE 0102.2— 7219. We retrieved the data from the 
Data Archives and Transmission System^ of ISAS/JAXA. 
All data were acquired with the normal full window mode 
and the 3x3 or 5x5 editing mode (Koyama et al. 2007). 
We used the data with the ASCA grades of 0, 2, 3, 4, 
and 6. As is mentioned in Koyama et al. (2007), a small 
fraction of the charge in a pixel is left behind (trailed) 
to the next pixel in the same column during the transfer. 
All data used in this paper were corrected for the trail 
phenomenon. The observations are summarized in table 2. 

Onboard Calibration Source Fe 

The calibration source ^^Fe produces the Mnl Ka line. 
The theoretical line center energy is 5895 eV (Bearden & 
Burr 1967; Krause & Oliver 1979). We used the data from 
August 2005 to April 2007. 

The Perseus Cluster of Galaxies 

This is one of the X-ray brightest clusters of galaxies 
in the sky. The X-ray spectrum is that of a thin thermal 
plasma with the strong Ka line of FeXXV. The plasma 
temperature changes smoothly from kT ~ 4 keV to ^ 
7 keV toward the outer region (Churazov et al. 2003), 
and the center energy of the Fe XXV Ka triplet is almost 
constant (~6.56 keV at z = 0.0176) within this tempera- 
ture range. Its radius of ~ 15' can cover the entire field of 
view of the XIS (18' x 18'). Thus this source is suitable for 
measuring the positional dependence of the energy scale. 

iS0102.2-7219 

This is one of the brightest supernova remnants in the 
Small Magellanic Cloud. With the spatial resolution of 
Suzaku, it can be regarded as a point source. There 
are many bright emission lines originated from thermal 
plasma in the X-ray spectrum below 2 keV. These lines 
are resolved with the XMM-Newton RGS, and the accu- 
rate energies of the line centroids are known (Rasmussen 
et al. 2001). This object has been used by many instru- 
ments for the calibration in the low-energy band, and an 
empirical model to describe the spectrum has been estab- 
lished^. We used this source as the energy-scale calibrator 
in the low-energy band. 

http://darts.isas.jaxa.jp/ 
^ http://cxc. harvard .edu/acis/E0102/ 
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4-2. Positional Dependence of the Energy Scale 

For the data of the Perseus cluster of galaxies, we di- 
vided the imaging area into four regions along the ActY 
axis, and extracted a spectrum from each region. Then we 
fitted the spectra in the 5-7.3 kcV band with a power-law 
model and a Gaussian function, and obtained the center 
pulse height of the FeXXV Ka line. Figure 5 shows the 
center pulse height as a function of i. Triangles and cir- 
cles indicate the data before and after the CTI correction, 
respectively. We can see no significant i dependence after 
the CTI correction, and this supports the validity of our 
correction. 

4-3. Energy Scale 

The goal is to determine a relation of PHo and X-ray en- 
ergy E. From the ground experiments, we found that the 
PHo-E relation can be expressed as a broken-linear func- 
tion linked at the Si-K edge energy of 1839 eV (Koyama 
et al. 2007). We then determined the PHq-E relation of 
each segment by using the lines of the calibration sources 
(Mm Ka line at 5895 eV) and 1E0102.2 - 7219 (Ka lines 
of O VIII , NelX and NeX around 650-1020 eV). 

We show the results after the CTI correction and the 
PHo-E conversion. Figure 6 shows measured center en- 
ergies of the Mnl Ka line as a function of time. Each 
mark in the plot has an effective exposure of more than 
60 ks. The mean values of the center energy are 5896.2, 
5895.4, 5895.0, and 5895.4 eV for XIS 0, 1, 2, and 3, 
respectively. The deviation around the theoretical center 
energy (5895 eV) is 7.8, 4.4, 6.6, and 7.8 eV for XIS 0, 1, 
2, and 3, respectively. Therefore, the time-averaged un- 
certainty of the absolute energy is ^ ± 0.1 % for the Mnl 
Ka line of the calibration sources. 

We also studied the time evolution of the deviation 
around the theoretical center energy, and the results are 
shown in figure 7. We can see that the deviation gradually 
increases with time. 

Figure 8 shows the center energy of the O VIII Ka line 
from the IE 0102.2 - 7219 data. The mean values of the 
center energy are 652.6, 653.8, 652.7, and 652.8 eV for 
XIS 0, 1; 2. and 3, respectively. The deviation around 
the center energy of the empirical model (653 eV) is 1.4, 
1.4, 2.3, and 1.1 eV for XIS 0, 1, 2, and 3. Therefore, 
the uncertainty of the absolute energy is ^ ± 0.2% for the 
OVIII Ka line of 1E0102.2 ~ 7219. 

4.4- Energy Resolution 

We examined the energy resolution in FWHM (A£') 
(Koyama et al. 2007) for each sensor; AE is common to 
all segments. We expressed AE as 

where a and b are time dependent parameters and AEq 
is the energy resolution determined by the ground exper- 
iments and obtained using equation 1 in Koyama et al. 
(2007). We determined a and b by using the time his- 



tory of the calibration sources and IE 0102.2 — 7219. The 
AE values obtained in this way is incorporated into the 
redistribution matrix file (RMF). 

Figure 9 shows the energy resolution of the Mnl Ka 
line after the column-to-column CTI correction. We also 
plot the results of the CTI correction, where we used the 
CTI values averaged over a segment (the column-averaged 
CTI correction). We can see that the energy resolution 
is greatly improved by the column-to-column CTI cor- 
rection. For example, the energy resolution in December 
2006 was greatly improved from ^^230 eV to ^ 190 eV. On 
the other hand, with the column-averaged CTI correction, 
the energy resolution is ^^230 eV and is not significantly 
improved. In figure 10, we compared the energy resolution 
of the Mnl Ka line with our RMF model. The deviation 
of the data points around our model is 5.6, 4.9, 3.4, and 
6.3 eV for XIS 0, 1, 2, and 3, respectively. 

5. Summary 

We have conducted the CFCI experiments six times 
in orbit. The CTI correction has been done with the 
CFCI results. We calibrated the energy scale of the 
XIS precisely using the onboard calibration sources and 
IE 0102.2 — 7219. Our calibration results have been ap- 
plied to all the data obtained with the normal full win- 
dow mode without the spaced-row charge injection. The 
results of the CFCI experiments and the current calibra- 
tion status are summarized as follows: 

1. We determined the CTIl and CTI 2 values of each 
column precisely based on the data of the CFCI ex- 
periments. We also found that the pulse height de- 
pendence of the CTI does not change with time. 

2. After the column-to-column CTI correction, we de- 
termined the PHq-E relation. We also modeled the 
time-dependent energy resolution. 

3. The uncertainty of the energy scale is ± 0.2 % for 
the OVIII Ka line (- 0.65 keV) of IE 0102.2 - 7219, 
and ± 0.1 % for the Mnl Ka hue (~5.9 keV) of the 
calibration sources. 

4. With the column-to-column CTI correction, the en- 
ergy resolution at 5.9 keV in December 2006 was 
greatly improved from ~230 eV to ~ 190 eV. 
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Appendix. Determination of the CTI values be- 
fore the first CFCI experiment 

First, wc determined the CTI values of the segments A 
and D in August 2005. We combined the data of the cali- 
bration sources from August 11 to 31, 2005, and obtained 
the PH' of the Mnl Ka hue. We also estimated PHo 
at 5895 eV from the PHq-E relation determined by the 
ground experiments, and obtained the ratio PH' / PHq. 
From equation 2, the ratio can be expressed approxi- 
mately as PH' / PHo = 1 — icaici — 1024c2, where icai is 
the mean transfer number of the calibration events (typ- 
ically ~900). We determined ci and C2 at 5895 eV from 
PH'/PHo with the Ci/c2 ratio fixed to the values shown 
in section 3. The ci and C2 for other pulse-height val- 
ues were calculated from equation 3, where we used the 
column-averaged and time-averaged /? values determined 
in section 3. Then for segments B and C, we took the av- 
erage CTI values of the segments A and D. We regard the 
CTI values obtained in this procedure as those on August 
11, 2005 (the day of the XIS first light). Note that these 
values are determined for each segment, not for each col- 
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Fig. 1. Time evolution of /3 averaged over one sensor. We 
sliow tlie standard deviation of /3 as the error bar. 
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Fig. 2. The number distribution of the CTIl parameter (ci ) . 
The parameter ci is obtained by the checker flag charge in- 
jection (CFCI) experiment on July 17, 2006. We show the 
result of XIS as a typical example. The typical error of ci is 
~ 4x10"'^. The mean and the peak (the most probable) val- 
ues are 9.4xl0~^ and 6.8xl0~^, respectively. The equivalent 
X-ray energy of the reference charge packet is 4.2 keV. 
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Table 1. Log of the checker flag charge injection experiments in orbit. 

Log Number Date Exposure The equivalent X-ray energy of 



(ks) the injected charge packets (keV) 











XIS 


XIS 1 


XIS 2 


XIS 3 


1 


2006/1/17- 


-20 


3.7 




7.1 




4.6 


2 


2006/5/20- 


-21 


1.2 


4.2 




3.9 




3 


2006/6/26- 


-27 


4.9 




0.3/7.1 




0.5/4.5 


4 


2006/7/17 




5.7 


0.6/4.2/8.0 




0.6/3.9/7.8 




5 


2006/8/25- 


-26 


4.9 




0.3/7.1 




0.5/7.0 


6 


2007/9/28 




25.3 


0.6/4.1/7.9 


0.4/7.2/7.9 




0.5/4.5/6.5 



Table 2. Observational log of the celestial objects. 



Object 


Obs.ID 


Date 


Exposure (ks) 


The Perseus cluster of galaxies 


800010010 


2006/2/1-2 


50.4 


1E0102.2-7219 


100014010 


2005/8/31 


24.3 




100044010 


2005/12/16-19 


59.7 




101005010 


2006/4/16 


21.3 




101005040 


2006/7/17 


20.6 




101005070 


2006/10/21-22 


18.5 




101005100 


2007/1/15 


22.6 




101005120 


2007/3/18-19 


18.2 




Ci at day 810/ Ci at day 372 

Fig. 4. Time history of column-averaged ci at 7.0 keV. We 
obtained ci at day 32 (August 11, 2005) from the calibration 
sources and that at the other day from the CFCI experiments. 
We used the relation cj oc (PHo)~^ when obtaining ci at 

Fig. 3. Same as flgure 2, but the ratio of ci obtained on 7.0 keV. Here we used (3 shown in section 3. 

July 17, 2006 (372 days after the launch) and on September 

28, 2007 (810 days after the launch). The typical error of the 

ratio is ~ 0.1. The mean and the peak values are 2.89 and 

2.63, respectively. Time evolution of Ci (increase with time) 

is clearly seen. 



Energy-Scale Calibration of the XIS 
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Fig. 5. Transfer number dependence of the center pulse 
height for the Fe XXV Ka line of the Perseus cluster of galax- 
ies. The result of XIS is shown as a typical example. 
Triangles and circles represent the data before and after the 
CTI correction, respectively. 



Fig. 7. The deviation of the center energy of the Mm Ka 
line from the theoretical value (5895 eV) as a function of time. 
Each mark (except for the last one of XIS 2) is obtained from 
forty data points of figure 6. 
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Fig. 6. Center energy of the Mn I Ka line for XIS 0-3 after 
the CTI correction. The theoretical center energy (5895 eV) 
is shown with dotted lines. Each mark in the plot has an 
effective exposure of more than 60 ks. 



Fig. 8. Center energy of the Ovill Ka line from 
1E0102.2 - 7219 after the CTI correction. Dotted lines in- 
dicate the center energy of the empirical model (653 eV). 
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Fig. 10. Energy resolution of the Mm Ka line for XIS 0-3 
after the CTI correetion. Eaeh mark in the plot has an ex- 
posure time of more than 60 ks. The solid lines indicate our 
model of the energy resolution which is incorporated into the 
redistribution matrix file. 



Fig. 9. Comparison of the column-to-column CTI correction 
(circles) and the column-averaged CTI correction (triangles). 
We show the energy resolution of the Mn I Ko line of XIS 
(a) and XIS 1 (b). Arrows in the plot indicate the days we 
performed the CFCI experiments. 



